Statistical Arbitrage and Optimal Trading with Transaction 

Costs in Futures Markets 



Theodoros Tsagaris* 
December 2, 2008 

Abstract 

We consider the Brownian market model and the problem of expected utility max- 
imization of terminal wealth. We, specifically, examine the problem of maximizing the 
utility of terminal wealth under the presence of transaction costs of a fund/agent invest- 
ing in futures markets. We offer some preliminary remarks about statistical arbitrage 
strategies and we set the framework for futures markets, and introduce concepts such as 
margin, gearing and slippage. The setting is of discrete time, and the price evolution of the 
futures prices is modelled as discrete random sequence involving Ito's sums. We assume 
the drift and the Brownian motion driving the return process are non-observable and the 
transaction costs are represented by the bid-ask spread. We provide explicit solution to 
the optimal portfolio process, and we offer an example using logarithmic utility. 

Keywords: Portfolio process, convex analysis, transaction cost, martingale representation, 
duality, utility maximization, statistical arbitrage. 

1 Introduction 

1.1 Optimal Investment 

This section examines the problem of a fund/agent with assets under management x, who 
wishes to invest in d assets, and aims to maximize the expectation of the utility function of 
the terminal wealth Xtj^ . The agent aims to invest her initial capital on d-dimensional futures 
contracts, which are driven by a d-dimensional Brownian motion. The manager /investor has 
only partial information, in the sense that she neither observes the drift /3, not the Brownian 
motion W. In addition, the investor incurs a transaction cost equal to half the bid/ask 
spread c for each lot(contract) traded. The investor's objective is to achieve an optimal 
position {Pt„) process that can hedge against the risk that is associated with the market, 
and deliver a constant level variance of the portfolio returns taking into consideration in the 
wealth equation any capital that is not required by margin calls. Finally, the manager would 
like to have the option to invest in a risk neutral manner i.e. the correlation between the 
underlying asset prices and the investment returns is 0. 
The following text is structured as follows: 
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We proceed in ll.3l with some explanatory remarks on relative/statistical arbitrage trading 
strategies commonly employed by proprietary desks. We then continue 12.11 on the definition 
of the futures price evolution process and the conditions under a unique and sufficient solution 
exists. We operate in the usual complete probability space {^,J^, {T)t^^Q, P), and our pro- 
cesses are adapted to the augmented filtration J^t„ generated by non-decreasing a-algebras. 
We define in addition the filtration generated by the futures prices and geometric returns J-f , 
respectively and show by direct observation that are equal. 

In section 12.21 we offer a remark of how to construct a statistical arbitrage strategy using 
the innovations approach, a tool very commonly used in stochastic filtering theory (see |2^)). 
Then, we continue with l2.3l where we offer necessary definitions and remarks on the operation 
of futures funds. We introduce some essential notions for the derivation of the optimal trading 
strategy such as gearing, position, unit value, contract size, margin and slippage. In [231 we 
define the relative risk process under transaction costs, and determine a new measure under 
which the expected value of the discounted final wealth is equal to the initial assets under 
management. We introduce a new formulation for the wealth process, and show explicitly 
from first principles how the wealth stochastic difference equation can be constructed using 
actual (bid/ask spread) rather than proportional transaction costs to the wealth process. We 
obtain an explicit solution for the cost term (see eq. ()19ll ) 
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and gain the insight that the cost term tends to infinity as the partition 11 = max^ At„. tends 
to 0. We also introduce the adjusted discount process as it has been amended for margin 
calls. 

For the remainder of the paper, we make use of the auxiliary results, presented in appendix 
[At regarding concave functional evaluated in some random process realization, and we state 
some known results about utility functions. We continue with section [3l where we establish 
the primal optimization problem and its dual. The main result is presented in section [U 
where we prove the optimality of the portfolio process vr in a similar manner to (15), and 
we characterize the optimal portfolio process using the martingale representation theorem. 
Finally, in section [5l we offer an example using logarithmic utility, and show explicitly that the 
terminal wealth maximizes the value function. Without sacrificing generality, and in order to 
offer an explicit solution to the reader, we suggest an approximating optimal portfolio process 
by using an approximate functional c for the cost term c. A remark is finally made that the 
investor should be interested in solution that satisfy the feasibility condition |/?| — |c| > 0. 

1.2 Previous Work 



The theory of optimal investment dates back to the seminal paper of Markowitz (l25l ) whose 



idea based on the examination of effective diversification of an investor's monetary resources. 



In the seminal papers of Merton (|27[ : [26l). a solution to the problem of optimal investment is 
presented for complete markets in the continuous case using tools of optimal stochastic control. 
Merton uses Hamilton-Jacobi-Bellman equation to derive a partial differential equation for 
the special case where asset prices are Markovian. 

In 1973, Bismut H) examines the optimal investment problem using a conjugate function 
in an optimal stochastic control setting. Bismut contribution forms what is known now as 
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the martingale approach in convex duahty theory. Harrison and Phska (|l2l ) provide a generic 
market model in continuous market trading and show that the optimal control problem can 
be solved using the martingale representation characterization of the terminal wealth by 
investing in complete markets. Cox, Huang (j^) present sufficient conditions for the existence 
of solution to the consumption-portfolio problem under a certain class of utility functions. 
Karatzas, Lehoczky and Shreve (fisl ) split the optimal investment problem into two more 
elementary problems where a small investor (one who cannot affect the market through his 
transactions) attempts to maximize the utility of consumption and terminal wealth separately. 

Further studies by Cvitanic, Karatzas p) generalize the results to the constrained case 
where the portfolio takes values in a predetermined closed, convex subset. Cuoco (0) extends 
the latter problem to the case of stochastic initial endowment. Karatzas, Lehoczky, Shreve 
and Xu (|l6|) study the case where the number of stocks is strictly smaller than the dimension 
of Brownian motions (incomplete markets) and suggest a solution by introducing additional 
fictitious stocks to complete the market. Xu, Shreve (jssl ) consider the finite horizon problem 
of optimal consumption/portfolio under short selling prohibition. 

Kramkov, Schachermayer (jiol ) study the problem of optimal investment in the more gen- 
eral framework of semimartingales, and provide necessary and sufficient conditions for the 
solution to the maximization of the expected utility of terminal wealth. They introduce 
the asymptotic elasticity of the utility function and show that it is necessary and sufficient 
condition for the optimal solution if it is strictly less than 1. 

Transaction costs are an important financial feature whose omission overestimates the 
wealth process, and subsequently affects the theoretical or derivative prices of the underlying 
assets. Davis, Norman (jlOj) have showed that the optimal policies are local times at the 
boundary, and that the boundary is determined by a nonlinear free boundary problem. Shreve, 
Soner (j34l ) are influenced by Davis, Norman work to provide an analysis based on viscosity 
solutions to Hamilton-Jacobi-Bellman equation and additionally show that for a power utility 
c^/p, < p < 1 there is a correspondence between the value function and transaction cost 
in the order of 2/3 power. Zariphopoulou OSg) shows that the value function is a viscosity 
solution, in the sense of Crandall, Lions (14), of a system of variational inequalities with 
gradient constraints. Cvitanic, Karatzas (|7|) provide an analysis of hedging a claim under 
transaction costs using duality theory. Cvitanic, Pham, Touzi apply to show that the 
cheapest strategy to hedging is a buy and hold strategy of a call option. Kabanov, Last (jl^ ) 
deploy the more general semimartingale model with transaction costs in a continuous setting. 
Under assumption of properness of the solvency cone, they provide characterization of the 
initial endowments that allow investors to hedge contingent claims. 

Lakner (23) approaches the optimal investment problem from the partial information 
case. He deploys some classical theory (Kalman-Bucy) to find an explicit solution for the 
portfolio process when the drift process is the unobservable component. Pham, Quenez (fsol ) 
use stochastic duality and filtering theory to characterize the optimal portfolio process subject 
to the assumption the stock prices follow a stochastic volatility model. 

1.3 Market Neutral Statistical Arbitrage Trading Strategies 

Statistical arbitrage methodologies refer to investment strategies that have high probability of 
outperforming the market (see notes in chapter 4 of (jH)). A very common, and particularly 
attractive to investors, statistical arbitrage strategy is called market-neutral strategy. Funds 
that deploy market-neutral strategies (sometimes refer to stat. arb. funds) typically pledge 
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to investors to deliver portfolio returns that are uncorrelated on average with the underlying 
asset or sector in question. 

Market neutral strategies are generally multivariate techniques (or bivariate in the special 
case of pairs trading) to define estimated equilibrium levels for a basket of assets, one may 
arbitrarily deploy. Such techniques generate a synthetic asset based on certain variables that 
could be interpreted as the fair price for a subset of factors from the available universe. Com- 
monly used variables are other assets classified in the same sector /industry or macroeconomic 
factors. Deviations from the fair price are considered as mis-pricings from the equilibrium 
price and one could use an arbitrary trading strategy to exploit such seemingly opportunity 
patterns. The forecasting ability of these strategies lies in the mean-reverting dynamics of the 
innovation process (for statistical arbitrage related remarks/techniques and an application, 
see B). 

In this study, we motivate our methodology on market neutral strategies and we offer 
some explanatory remarks of how such strategies can be constructed using standard filtering 
techniques. These techniques provide estimates of drift that are uncorrelated to the underlying 
price, which in turn are deployed in the optimal portfolio characterizations derived by classical 
duality theory in optimal investment. 



2 The model 

2.1 The Market Model 

Assume we are given a complete probability space {Q, T, {T)t,^>o, P) , and let J^tn be the 
augmented filtration generated by the union of a non-decreasing family of sub-a-algebras and 
AA-null sets. 

Let T be a finite closed interval [0, T] , where T is the terminal time. Assume we have the 
following partition = to < ti < . . . < tn ■ ■ ■ < tj\f = T oi the set T, and we denote the set of 
non-decreasing elements by 5 = {tQ,ti, . . . ,tn, ■ ■ ■ ,t]\f}- Moreover, we introduce the following 
notation for the time increment At = At„ = tn+i — tn, and assume equidistant increments 
tn = to ~^ ^n, where 11 is the maximum timestep 11 = max^ At„. 

We introduce the process Ft„ as a mapping from J7 x 5 — > M*^, for which is measurable with 
respect to the smallest-sigma field of subsets oi Q x S. In this setting, there is no meaning 
to assume right continuity. The optional c-field O is the fi-field generated by all adapted 
processes to {Tt„)neN, such that Ft„ is Tt„ measurable for each n G N. 

Denote by 

J^f^=a{Ft„;tneS,neN) (1) 

the fj-field generated by the futures contracts prices. The sigma field is a subfield of 
{^t„]in ^ S,n G N}. The filtration represents the information that is available to 
investors up to, and including time 

We adopt a financial market consisting of d-futures F^'^\uj,tn)-, ^ = 1, Vt„ G S that 
satisfies the following evolution equatiorE] 

Ai^t? = ^*? f /3«(UAt + X:a(^^')A<-M (2) 



^Henceforth, we also make use of the following condensed notation: /t„ — f{Lu,tn). 
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where 



• Wt„ is ad dimensional Brownian motion with correlation matrix p i.e. £^[APFt^ AW^* ] = 
pAt. 

p = p^^^\ i = 1, . . . ,d,j = 1, . . . ,d is a d X d non-random matrix satisfying p^'''^ < oo 
(where (•) denotes summation over all subscripts) 

• Pt„ is a d dimensional mean rate of return process subject to Ylk=o Wf^tkW^^ < 

Pt„ : ^ X S R'^ is a measurable mapping such that Pt„ is J^tn measurable for every 

• a = a^^^\ i = 1, . . . ,d,j = 1, . . . ,d is a d x d non-random volatility matrix satisfying 
cj^'''-' < oo 

• The initial prices is a d-dimensional vector of deterministic constants 
and the following statements hold 



and J^t„ is independent of c7(AWf„) 



Assumption 1. We assume a and p have full rank, therefore are nonsingular and hence 
invertible. (See Lemma 2.1 of Karatzas, Lehoczky, Shreve (1987) on uniform ellipticity as 
sufficient condition for invertibility ) 



By Karatzas and Shreve, eq. ([2]) has a weak solution (see (llTI ). Chapter 5, Problem 6.15). 
All vectors are column vectors, unless stated otherwise. We will denote the transpose by 
* (i.e. (v*)* = v). 

We assume that the agent can only observe the futures prices but the drift process and 
the Brownian motion W are unobservable. The initial prices F*''^^ and the volatility matrix a 
are known in advance. 

Moreover, we denote the quadratic variation of Wt„ by 



n-l 

[w,wu = - Wt,){Wt,^, - Wt, 



k=0 

Assumption 2. We will be making use, throughout the remainder of the text, of the simpli- 
fying assumption that the the quadratic variation can be replaced by its expectation. 

Finally, throughout the text, we assume that any sequence of measurable mappings 
Yn, n G N, takes values in a Polish space {E,B) and hence we have that for any Borel 
function /i : S ^ M , the function ^(1^) converges to h{X) surely. The above assumption is 
crucial for the existence of conditional probability of a mapping of Y given a sub-cr-algebra 
ECB (see Stroock,Varadhan 2000). 

2.2 The Trading Strategy 

Define the d-dimensional return process R : S M.'^ by 

^p{i) ^ p(*)Ai?« (3) 

Lji iJir i^n ^ ' 
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and applying ([2]) in the above, we obtain 

d 

AR^^ = /3f ^ At + y a^'^^ Awi'^ (4) 

3=1 

Lemma 3. The filtration generated by the return process 

= a{Rt, Rt„ } - null sets of ^} 

is equal to the augmented filtration defined in ([7P = 

Proof. Eq. ^ implies that both TJ^, J-^ generate the same filtration. □ 

From the above lemma and for the remainder of the text, we take the return process 
R = {Rt„,TiD as the observable random process, and /3 = {Ptn,^t„) as the unobservable 
component. 

Assumption 4. The linear dynamics of f3t„ are given by the d-dimensional process 

d d 

j=i j=i 

where W^'^^ is a d-dimensional Brownian motion and W , T^(^) are independent. The matrices 
a and q are non-random in R'^^'^ and satisfy the following conditions 

n—l n—1 
fc=0 k=0 

Pto is the initial condition satisfying i?[||/?to || At] < oo and is independent of cr{Wti , ■ ■ ■ , Wt„ } 
and C7{<') 

(2) 

Moreover, Tt„ is independent of a{AWf:^ ). 

This is a classical partial information estimation problem of finding the state process f3t„ 
by observing the returns of the futures prices. The investor only observes the futures contract 
prices, not the Brownian motion or the drift. For a broad class of mean square cost functions, 
it can be shown that the optimal estimate, in terms of minimising the expected value of the 
cost functions, is the expectation conditional on the filtration generated by the observation 
process. 

Let pj^ denote the projection from the Hilbert space L'^{P) onto the closure in Lp'{P) of 
linearly combined functions 

N = {ao + aiRj^ + ... + adPf^ G -^^(-P) ; -R is measurable} 
Then, coincides with the conditional expectation 

P}^^=E[Pi^\^Pj=M^ (6) 

(see (Q; 0; [El; S S 0) for proof of this known result). We denote by Pt„ the vector 
with elements {/d} ,/??,..., pfj. 
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Then, the innovation process is defined as 

Az.,„ = (a)-i(Ai?i„-A„At) (7) 

Az^t„ = iar\(3t,. - (3tJAt + AWt„ (8) 
Lemma 5. The conditional expectation Afj^ is given by 

and the conditional variance is of the order of 

Var[Ai.£VtU = Piiitn+i - tn) = At 
Proof. The first statement fohows form iterated conditional expectation 

n—1 n—1 

EWtjFt,]=E[^{a)-\p,,-A,)At + J2^W,,\Ft,] 

k=0 k=0 
n—1 n—1 

EWtjFt,] = EiY^i^r'Em, - paFf^WFtj + i^E ElAWt^lFf^Wt,] = 

k=0 k=0 

By discarding the quadratic At terms, and replacing the quadratic variation by its expecta- 
tion, we get 

n-l 

VarKJ^,lj = EiulutjTt-,] = E[E[Y^ A[^,;mj|Fa!F,lj 

k=0 
n-l 

VarK-z/tJ^,U = EpAt 

A:=0 

□ 

Remark 6. From lemma ([^, we deduce the mean reversion properties of the innovation 
process Av around a constant level. The stationary nature of the spread (innovation) has 
been widely exploited by proprietary traders to establish appropriate trading rules to benefit 
from the "mis-pricings" of the futures contract price. The futures contract price is commonly 
interpreted as expensive when the spread is positive and cheap when negative. 

The use of spread, for the establishment of proprietary trading strategies or hedging, is 
particularly attractive because is uncorrelated to the futures prices. By construction, the spread 
is uncorrelated to futures prices because of the orthogonality condition of the projection theo- 
rem. Any proprietary strategy the depends solely on the spread is typically a "market neutral" 
strategy because it is not influenced by the fluctuations in the underlying asset (market). 

In reality though, the proprietary choice of filter to calculate [3 leads to estimation issues: 
It introduces estimation error, as generally is notoriously difficult to capture the drift rate 
accurately, because of the low signal to noise ratio of financial data. And, moreover, the nec- 
essary and sufficient condition of orthogonality, may not always hold, and hence the strategy 
might be occasionally affected by market fluctuations. 
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2.3 The Investment Process 



In futures markets, an agent enters the market by paying usually a locally denominated 
amount that reflects the price of the contract at time t. Generally, the price of the contract 
is calculated by Ct^^ = diag(/)Ft^ ,where / E represents the unit value of the futures 
contract H and C : S ^ M^j. is a measurable process with respect to Tt„- The agent is 
required to set aside an initial and maintenance margin m € [0, 1] (henceforth to be taken 
as constant), generally 20% of the current market, unlike stocks where the margin could be 
roughly as much as 50%(sale proceeds plus 50%) for long (short) positions. 

Definition 7. Consider the mappings tt : S x n ^ R'^ and P"" : S X n ^ M'^. Then, the 
wealth process corresponding to the investment strategy vr and initial endowment x is denoted 
by X^^^ , and represents the income accumulated up to, and including time t. The portfolio 
or investment process vr represents the portfolio weight assigned to asset i (or correspondingly 
the position P^^ of futures contracts, as we will later establish). The portfolio process is 
measurable with respect to J-^, i.e. the investor's portfolio allocation is based on information 
received only by the futures prices, not the drift or the Brownian motion. 

In the following remark, we offer some explanatory comments for the portfolio process 
and the orderbook that would facilitate the build-up of the formulation in the following text. 

Remark 8. Consider at any point of time t, there is an orderbook of futures positions 

(ii) 

,j € 1, . . . J available, that consists of a stack of positions that reflects the bid and ask 
prices, where j refers to the ordering of the stack. We assume that the position, offered for 
purchase or bid for sale at the top of the stack 'pI^\ represents the narrowest spread c > 0. 

An investor whose preference is to take position PT^'^ > 'Pt^^ will need to move deeper in to 
the stack, and pay the additional slippage incurred^ 

Assumption 9. The position we take in the market is bounded from above P^^ < Vl^ , in 
the sense that we can not exceed the available position that is associated with the narrowest 
spread. 

Assumption 10. The assume that half the bid/ask spread c/2 represents the deviation in 
market unit terms between the intended and the actual "fill" price for any buy or sell trade. 
For simplicity, we assume that c is a constant and does not fluctuate during the period of the 
normal trading hours. 

The cost an investor incurs in the market is only equal to the slippage and there is no 
commission cost associated with any trade. This assumption is not restrictive as commission 
tends to be only a small contribution to the overall cost. 

Discussion 11. Contrary to the literature, where the optimal investment problem is set up 
on d + 1 assets, (e.g. d stocks and a bank account), it is reasonable for futures investment to 
consider that we invest the entire capital in risky assets, and that the proportion of capital 
1 — m, Mm G [0, 1] that is not required by margin calls is deposited at the risk-free rate. Note 
that the same capital is invested concurrently on risky assets and risk-free deposit account. 

^The locally denominated amount that corresponds to the one unit movement in the future market price. 
^For real time examples of the futures orderbook, check for instance Chicago Board Of Trade web site. 
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Consider an investor who wants to invest in d assets, and his initial endowment is x > 0. 
The portfoho process is an -measurable adapted process that satisfies 

k=0 

Note, that we allow for short or long position, and that we could relax any constraint on the 
position, subject to the fact that assumption ([9|) holds throughout this text. 

2.4 The Problem Formulation 

In this section, we use throughout for consistency the same notation as in (fl5l).(fl^). (flil ).(|6l). 

The next lemma is rather crucial for section continuity, as we establish the new measure 
under the underlying probability space. The new measure adjusts for the drift in (jl7p . and 
under this measure the wealth equation will turn into a local martingale. The transformation 
is attained through a representative of the Radon-Nikodym derivative 

dPn = ZtJPn (9) 

where Z is given in eq. (jlip , and assume it satisfies the Novikov condition (jl4p . Consequently, 
we establish the absolute continuity of measures Pn « Pn, and we deduce that 

n-l 



k=0 

is a Brownian motion under the new measure P„. 

Definition 12. Consider an investor who operates in the futures markets. The relative 
risk process under transaction costs is a progressively measurable process 9t„ : <S — > M*^ that 
satisfies 

P^^t„ Q„ (10) 

Definition 13. Let Zn he a Pn-martingale given by 

jt— 1 ^ n— 1 

Zt^ = exp{- O^AWis) - - ^ 9lpet^At} (11) 

fc=0 A:=0 

We denote the optional projection of Zn{t) to T^^^ by 

Ct^ = EiZtJJ^H] 

Lemma 14. Let 

n-1 

Wt„ = Wt„ + J2p0t,Atk (12) 



k=0 

where 0t„p9t„ satisfies the following condition 

n-1 



P[u;:^^r„p0i„Atfc<oo] = l 



k=0 
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If E[Zt„] = 1, then Zt„ is a local martingale and a supermartingale and there is another 
probability measure Pn that is J-t„-adapted, and corresponds to the Radon-Nikodym derivative 
dPn = Zt„dPn and Pn « Pn- The process Wt„ is a Brownian motion under the probability 
measure Pn- 

Proof. See appendix ([C]) □ 



E\m\\]<oo 



Theorem 15. Assume that 
Let 



Then, it follows that the stochastic difference equation 

A ( — ) = —e*, AWt 
holds, and its stochastic integral representation is given by 

n—l , n— 1 



= exp{- OlAWt, + IY1 ^Ipet.At} (13) 



2 

fc=0 k=0 

Proof. See appendix dC]) □ 
Definition 16. The financial market A4 is viable, if the following conditions are satisfied 

n-1 

P[uj:J2O*t,p0t,As<^] = l (14) 

A:=0 

Remark 17. The intuition behind the above argument is, that there is a price of risk that the 
market will quote you, that reflects the market risk net of transaction costs. In this setting of 
futures investments, and contrary to the literature, it is economically intuitive to formulate 
the market price of risk without the usual risk-free rate term that translates the market risk to 
excess market risk. As suggested in remark (li ij). in futures trading, the investor earns proceeds 
from the capital twice, through investing both on risky and risk-less assets concurrently. On 
the other hand, a long-short strategy in stocks is cash-intensive, and the same capital has to be 
either allocated in the risk-free account or on some risky strategy, which makes the trade-off 
between risky and risk-free assets relevant. 

Define the adjusted discount process by 



n-1 



7i„ = exp 
and the state price density by 



^{l-m)rAt\ (15) 



k=0 



Ht„=lt„Zt„ (16) 

and recall that m (see remark (jlip ) is a non- negative constant with compact support in [0, 1] 
that represents the margin as proportion of the wealth X^^^ . We also introduce the following 
notation 

^tn = 7t„Ct„ 
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Definition 18. The initial endowment x > is independent of the a-algebra generated by 
F{-), and satisfies 

E[\\xf] < cx) 

A portfolio process tt{-) is called admissible if it satisfies 

X > and X^'"" > Vt G 5 
The set of all portfolio processes, where the above conditions hold, will be denoted by A{x) 

Assume there is an agent, who can decide on a non-anticipative basis, to invest her capital 
at time t„, to the d available futures by selecting an appropriate vector of positions P^^. We 
assume that the positions are sufficiently small, in the sense that they do not "move" the 
market (see assumption [9]) . Then, the wealth equation of an investor X^^^ with portfolio 
process vfj^ and initial endowment x, is derived as follows: 

AX^f = (1 - m)rX^fAt + AFl di^g{f)P:f -\e diag(/) | AP,');'^^ | 

AXl^^ = (1 - myXl^^At + Xl^^AFl diag(/) diag(C,;iX - \e diag(/) | AP,:'^^ | 
AXf:^ = (1 - m)rXr:^At + Xf^^AP,; diag(P,;i)<„ - \e diag(/)| AP.'^'^^ | 

AXf;^ = (1 - m)rXf;^At + - \c {{P^fT diag(C* JAt)"^ diag(/)|AP,:'!^ \)At 

+X-/«)VAt^,„ 
AX^f = (1 - m)rXl^At + X^fininPt,. - ctJAt 

+ X:f{niraAWt^ (17) 

where 

P-^ = X,ydiag(Cr;)^t (18) 

is the futures contacts position. The "one-way" spread relative to the price is a c : SxQ 
process and is given by 

c«/W|AP;'^(*)| 

2p;''=«d*)At 

for the ith asset at time t. We denote the vector formulation of c by q„ = {cj^, . . . ,cf^}, 
and we refer to | • | as the component-wise absolute value. The superscript in ir^ denotes the 
dependence of the portfolio process to non-zero transaction cost. 

Remark 19. Let k be the gearing /leverage weights that take values in M^. In managed 
futures funds, the gearing vector is a proprietary vector of weights that plays the important 
role of establishing the leverage of the fund. A typical futures fund pledges to investors to 
fluctuate within a certain band of volatility of monetary returns^ Generally, futures funds 
tend to market their strategies between 10%- 25% of annualised volatility. 



''We define monetary returns in the interval in], ft G N by ^J;^^^^ 'YlZ^n-hi^t^'^y diag(/)A_Ft^ 
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Discussion 20. The first term in eg. JiTp refers to the capital earned by the wealth invested 
in the interest free account. Note, that we adjust the wealth by 1 — m factor, as we assume 
that we do not earn interest on the margin. The second term refers to capital invested in the 
d-dimensional futures contracts, where tt refers to the portfolio process. Finally, the last term 
is the diffusion term of the wealth equation. Note that the position P'"'^ can be arbitrarily 
multiplied by the leverage term k in order to achieve the required volatility of returns a 

Var ]=a, aen+ 

\ ^tn~At J 

where At typically represents a daily time interval. 

From eq. (|12p and l\W \ 1171 [T2]) . we deduce that the wealth equation permits the following 
stochastic representation 

AX^f = (1 - m)rX:fAt + Xf/«J*(A„ - ctJAt 

+X^f{7rira{AWt„-pet„At) 
= (1 - myX^fAt + X^finlTaAWt^ (20) 
Lemma 21. The margin adjusted, state price density portfolio process 

-Yt Zf X, 

I t-n tfi 

is a local martingale and a supermartingale. 

Proof. Note that 7f„ is a decreasing function of t. From the monotonicity of jt„ and the 
Novikov condition (jl4p . it follows that jtnZt^X^^^ is bounded. 

Then, by direct application of Taylor theorem and discarding the higher order terms of 
At, we obtain 

A^t^X^f = -7i„(l - m)rX^fAt + 7t„(l - m^X^fAt 

^,„X^f = x + J2lt,X^f{7r^*Am, (21) 

k=0 

Now, recall that H^^^ is the state price density given by (|16|) . and apply Taylor expansion 
in the AHf The boundeness of AHf XT'^' follows using similar arguments as above. 

Arn^X^f = (A7i„Xf;^)Z,„ + (AZ,j7,„Xr;^ + A[7^„Xf;^ Z^J (22) 

Note that AZt„ = —9t,^Zt„AWt„ follows from Taylor's theorem application (we already 
established that in lemma (fl^ 



Ht„X:f{7TiraAWt^ - Ht,,X^fdt^AWt„ 

-Ht^X^finlYapet^At 
= Ht„Xlf{a7rl-etSAWt^ 
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or equivalently in Ito's sum form 

n-l 

m^x^f = x + Y^ m^x^f {aTTi - e,^y aw,, (23) 

fc=0 

Eq. (j23p is a stochastic integral, and therefore a local martingale. Subject to the Novikov's 
condition, and by Fatou's lemma, we deduce that any non-negative local martingale is a 
supermartingle. 

It follows from above that the budget constraint equation is given by 

ElHt^X^f] < X (24) 

□ 

To ensure text continuity, we refer the auxiliary results about utility functions stated in 
the section IA.2I of the appendix, and we introduce a special class of normal concave utilities 
of stochastic processes, following similar arguments to that of normal convex integrands of 
Rockafellar (|3ll ) (see section lA.ll in the appendix) . 



3 The Optimization Problem 
3.1 The Primal Problem 

Consider an investor who wants to maximize the expected value of his terminal wealth 
E\U {X^^)] over the class of admissible trading stategies A (see def. p^). where U : (0, oo) — > 
M is a strictly increasing, strictly concave function, continuously differentiable utility function 
that satisfies 

ii;[C/(Xf/)] < oo (25) 

We name the objective function under the primal problem by P, and we denote the value 
function by 

V{x) = suvE[U{X-,;^^)] (26) 

subject to V{x) < oo, Vx € (0,00). The objective function ([26]) is deduced from (|15]) . when 
-in the notation of section dS])-, we evaluate U{y) at y = 0. The optimal value of P is 
therefore 

inf P = \B.iJV{x) 

The P can be less rigorously interpreted to the following dynamic optimization problem 

supi?[t/(xr/)] 

s.t. is financed 6y vr G A{x) & with the initial endowment being x 

Note that the above optimization problem has infinitely many budget constraints. 
A necessary and sufficient (for justification of the assumption, see KLSX fjld ) ) condition 
for the above assumption is 

U{x) <ki + k2x\ ki > 0,k2 > 0,5 e (0,1) 
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A portfolio process tt that achieves the supremum of the expected value of the utility 
function evaluated at the terminal wealth value at time T, is called optimal. 

Lemma 22. The value function V i26\) is a continuous, non- decreasing concave function. 

Proof. The non-decreasing property of V is easily observed by the fact that the set A{x) 
increases as x increases (see def. (jlSp ). Now, letO<p = l — (7<1 and xi,X2 > 0. From the 
concavity of the utility function U, we have 

pV{xi) + qV{x2) 
< V{pxi + qx2) 

i.e. the graph lies above any of its chords. Continuity is implied by concavity. 

□ 



3.2 The Dual Problem 

The dual portfolio process is an M'^ value process adapted with respect to that satisfies 

n-l 

E[J2 W^tkW^^t] < a.s. 

We name the objective function under the dual problem by D, and we denote the value 
function by 

V{y) = mlE[U{Y,Y)] 

and we assume that V{y) < oo, Vy G (0,oo), where A is the admissible class under the 
dual problem. The optimal value of D is 

sup D = sup V{y) 

Existence and uniqueness of the dual value function follows from the strict convexity of 
U{-) (see XSdH)). 

Lemma 23. The value function V ^26\) is a continuous, non-increasing convex function. 
Proof. The proof follows from similar arguments as lemma (j22p . 

□ 

Moreover, the following assumption (see (jiol )) ensures existence of the optimal solution. 



4 Optimal Portfolio Process 

In the following section we present our main result. We show that by solving the static 
optimization problem (see section |B] for a similar proof to (3) that the solution to the static 
problem coincides to the solution of the dynamic problem) we obtain an optimal portfolio 
process characterization. 
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4.1 Portfolio Process Characterization 

In the following section, we prove that the dual portfolio process yf provides a bound for the 
value function V. We follow similar arguments to KLSX ifl^ ). CK([^ to show the existence of 
the optimal 7r(-) under the primal problem. 

We introduce the following mapping X : (0, oo) (0, oo), and we define 



A:{y) = E[Ht^IiyHtJ], < y < oo (27) 

Note, that the monotonicity (strictly decreasing) and continuity properties of X{y) follow 
directly from properties of /(•). Moreover, define y{-) as the inverse of X(-). 

Theorem 24. Let 

^ = I{y{x)HtJ (28) 

and assume that satisfies the following condition E[U{(,))] < oo. The value function V attains 
its supremum if and only if 

mtA] = X 

Then, Vvr € A{x), we have 

V{x) < E[U{C)] 

Proof. See appendix (jC]) □ 



From theorem (j24p follows that there is a portfolio process vf"^ that attains the supremum of 

'TT,X /- 



the value function, and we call it optimal when we have the following equality X^'^ = ^tpj = 



for ^ defined in eq. (j28|) . 

The following lemma is in the spirit of (0; 15; 16: Issl). 

Lemma 25. The optimal portfolio process vf is defined by 

K = T^i'^yH9t„+MM (29) 
Proof. Define the positive martingale 

M{t) = E[h^xf;;'\J^tJ 

Then, there exists a unique stochastic process g{s,u;), such that Mt„ admits the following 
stochastic integral representation (see thrm 4.3.4 of Oksendal ijigl )) 

n-l 



Mt„=x + J29l^Wt,, a.syt>0 



k=0 



and the mapping g : S x n R'^, satisfies J]fc=d llftJP^* 

< oo. 

Firstly, note that from Taylor theorem and eq. (112^ . we have 



> /- I'll ''^ 



where 9t^ = {pa) ^(3t„- By Taylor expansion, we obtain 



A— Mi„ = A— Mi„ + (AMiJ— + A[— ,MtJ 
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= ^{gt„+MtJt^yAWt„ (30) 
Consequently, the discounted optimal wealth is defined as 

itA = ^ + Y.r (^*^ + ^tAy^^t>= (31) 

Substitute the optimal portfolio process (p9|) in ([2T]) to obtain the portfolio process given by 
(f3T]) . That completes the proof. 

□ 



5 Examples 

Consider the case, where the utility of the value function (|26|) is logarithmic U{x) = log{x). 
Then, according to theorem ()24p . the explicit formulation of the terminal wealth that maxi- 
mizes the value function is 

f n—l n— 1 ^ 

e = xexp ^^(1 - m)r + Jlpl^At + ^ OlAWt, (32) 



. fc=0 k=0 



where X{y) = ^ and y{x) = ^ 
Now, rearranging (|23p . 



X 



^ n—l 

^tf = ^ + ^ E H^^^tf - ^tj* AW^*, (33) 

and substituting the optimal portfolio process p4p 

= (^P^)"'(A„-Qj (34) 
Pff = Xff di^giCi-:)nl (35) 

into (p3]) . we conclude that 

Therefore, the optimal portfolio process is given by eq. (|34p and the value function is 

n-l 

y(x) = E[\ogX^;^] =\ogx + ^^(1 - m)r + ^^^p^t, Ai 

fc=0 
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Equivalently, when the bid-ask spread is zero, we have 



= {'yp'^r^Pt. (36) 

Pf„=xf;^diag(C,;i)?f,„ (37) 

which is simply the formulation of the optimal portfolio process/position without transaction 
costs. 

Remark 26. Substituting the optimal portfolio process \34^ into the position process 
and introducing a scalar multiplier, the gearing process (see remark ^19^) ). we have 

P!f = xffdiag{k)die.g{CrJ)nl 

Moreover, in order to retain the positivity of the expected payoff of the investor, we need to 
satisfy the following condition 

1^1 - |c| > 

or equivalently to introduce the following transformation 

Tt„ = max(/3t„ - q„ , 0) + min(^t„ - q„ , 0) (38) 

and hence, we obtain the following position and portfolio process 

Pff = Xff diag(fc) diag(C,;i)C (39) 

C = M-'rt„ (40) 
by using the following approximation for the cost term 



^i) ^ J \ tn tn-1 



(41) 



The final position functional ( eq. \39\) makes use of the approximating function c of c, where 
P^'^ is approximated by P^ , the optimal position process under logarithmic utility without 
transaction costs. 

The position process P^^^ is interpreted as the number of positions for each asset i subject 
to the gearing vector k, in order the portfolio to achieve the required volatility of monetary 
returns (see section i2.3\) for further explanation). The vector k could also be interpreted as 
a vector of proprietary portfolio weights. 

Discussion 27. The optimal portfolio/position process in ^39\) is a classic mean-variance 
criterion, with the difference that our P^^^^ mapping refers to the d-dimensional position we 
take in the futures markets. Moreover, contrary to the literature that reports ir as a propor- 
tion, in our context it represents an arbitrary weight, and moreover we note the introduction 
of few additional terms, namely the bid- ask spread, the contract size, the gearing, the correla- 
tion matrix, and the futures unit value(compare to (16) who derives the usual mean-variance 
criterion). 

The incorporation of transaction costs, has a remarkable effect in the position. We observe 
that the (3t„ term in l^34\ ) is adjusted by another process, namely ct,^ , that corresponds to the 
cost of entering that position. The latter process is simply the slippage (see def. 110\) ) incurred 
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in the market proportional to our expected movement. We explicitly assume a linear relation 
between slippage and magnitude of market order, but that is again reasonable subject to the 
assumption (^(See if7fl] for non-linear formulations). 

We need, additionally, to observe that at infitesimal time changes, our expectation will 
hardly change but the cost of purchasing /selling a contract is invariant of time changes (note 
that we only divide the cost term by At). 

Remark 28. Assume that the net return process Tt„ in i38\) represents our expectation of a 
particular trading strategy for asset i. Then, using Chebyshev inequality, we have that 

P{'Tt^>e)<^^, 6>0 
for significantly small e. But as 11 —> 0, we assert that Tf^ =0 i.e. 

P{^t„ >e)=0 

unless the increment At is sufficiently large such that \Pt„ \ ^ Q„- 



6 Conclusion 

We studied the problem of optimal investment under transaction cost and offered a number of 
results pertaining a relative arbitrage trading strategy. Our study focuses in futures markets, 
contrary to the literature where the bulk of the studies in optimal investment are in eq uities 
markets and few recent studies under transaction costs in currency markets (see (jil: [l3l)). 

The application of optimal investment strategies in futures markets necessitates the intro- 
duction of concepts such as position, gearing/leverage, slippage and margin, and the modelling 
of those in the wealth process. Moreover, we offer a brief description of the orderbook, and 
how it associates to the slippage used in the model under some boundness criteria. 

We show explicitly, through the construction of the partial information dynamics, how 
the filtration generated by the futures prices is equal to the filtration generated by the return 
process. We introduce the concept of relative arbitrage in systematic trading, and we offer 
some intuitive remarks how is related to the innovation process and why such strategy may 
be profitable. 

We proceed by constructing the wealth process by including transaction costs, and de- 
ploying the futures contracts position as control. These formulae account for the actual 
transaction cost, rather than deploying proportional constant transaction costs to the wealth 
process. Then, we use an arbitrary gearing factor to achieve the required rate of volatility 
the fund strategy entails. 

We introduce the normal concave utilities, in line with Rockafellar (jsil ). whose paper 
provides a detailed exposition of integrals which are convex functionals. Under this very 
generic class of utilities, we formulate the primal and dual problem, and provide conditions 
under which there is no duality gap. We continue using arg ^ments adapted from g 0; 0) 



to show that the value function attains its supremum under the static problem, which in turn 
proves to be equivalent to the dynamic. We conclude the section by providing an example 
using logarithmic utility, and an explicit representation of the optimal portfolio process. 

The above theory could be readily extended to the constrained case, in the spirit of the 
work of , as the position could be constrained due to liquidity issues. Moreover, the above 
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analysis of the trading strategy could encompass results in stochastic filtering theory to find 
estimators for the drift of the wealth process, and to study the properties of the estimator 
using statistical inference tools for stochastic processes. In addition, we could model the 
orderbook using stochastic processes and relax the assumption [9] that relates to remark [H] 
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APPENDIX 



A Auxiliary Results 

A.l Normal Concave Utilities 

In this section, we provide some results on concave utilities functionals evaluated at some 
random process realization for some elementary event cj, to serve as auxiliary results for 
establishing existence and uniqueness of solutions. 

Definition 29. A concave utility U{Xt), defined on [0, oo) and maps to M-|_,is normal if it is 
proper (i.e. epigraph is non-empty and contains no vertical lines) and is continuous, and if 
there is a family E of random variables X from [0, T] ^ M that satisfy: 

• VX^ G E, {Xt,J^t) is progressively measurable process i.e. for every time t, the map 
[0,T] xVL is B{\^,T])^ T measurable and adapted to Tt . 

• \/t £ [0, T], E f^doniU is dense in domU 

• € [0, T], U is essentially smooth, and satisfies the following conditions for C = 
int(dom U)) 

1. C is not empty 

2. U is differentiable everywhere in C 

Lemma 30. Let U be a continuous, proper concave functional in M+. Then, U is a normal 
concave utility. 

Proof. Since the effective domain of C/ is a nonempty concave set, we are able to define Z? as a 
dense subset of dom U . Moreover, let E be the family of random variables X from [0, T] — > R. 
Then, the conditions stated in def. (j29p are satisfied. □ 

Lemma 31. Suppose U is continuous in X for each t, and for each X is measurable in t, 
and has interior points in its effective domain. Then, U is a normal concave utility. 

Proof. Let the the family E of random variables take values in the dense subset D. Then, 
U is continuous in X , because D has dense intersection with domt/, and the measurability 
condition holds by hypothesis. □ 

Lemma 32. Let U the conjugate of the normal concave utility U. Then, for every progressively 
measurable process Xf from [0, T] M, the function U{Xt) is measurable. 

Proof. We have by definition that 

-U{y) = inf{< Xt, y > -U{Xt)\y € M} 

where the notation "<>" here denotes inner product. In the above formulation, we only need 
to show that any value of U {Xt) =< Xt,y > could be approximated by values in Ef] dom U. 
By hypothesis, we have that E P| dom U is dense and moreover given that the dom U is the 
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closure of the relative interior, the intersection of E with the relative interior must be dense. 



According to (j32l ). U is continuous in the relative interior of the effective domain U, and the 
values at boundary points can be obtained as limits. Hence, the values of U are limits of 
those for E P| dom U. 

It follows that the pointwise infimum is a collection of measurable functions, and therefore 
measurable. 

□ 

A. 2 Utility Functions 

A utility function U is a strictly increasing, normal concave mapping that satisfies the Inada 
conditions 

[/'(0+) = lim U'{x) = oo, U'{oo) = lim U'{x) = (42) 
and the following growth condition 

<U{x) <a{l + x'') (43) 

for a > and ly G (0, 1). 

Denote by / : (0, oo) (0, cxd) the strictly decreasing, continuous inverse of U'. This 
function satisfies 

7(0+) = oo, /(oo) = 

and 

U'{I{y)) = y Vy > 0, U'{I{x)) = x Vx > 

Let T and T be two Hausdorff topological vector spaces separated by a closed affine 
hyper plane Ti, which can be considered as bilinear form over T y.T . Let C C T be a convex, 
closed and non-emptv set and denote by D the dual of C (see Hahn-Banach Theorem in 
Ekeland, Temam (I 111 ) for existence of non-zero continuous linear forms in T). This is to 
say that any continuous function of x can be expressed as a linear function of < x, y > and 
conversely. 

Lemma 33. Let U he a normal concave utility. Then the Legendre conjugate is defined by 

U{y) = U{I{y)) - I{y)y (44) 

Hence, the Legendre conjugate (D, U) of (C, U) is well defined, where D is the image of 
C under the gradient mapping \/U. 

Proof The sub differential mapping 'dU reduces to the gradient mapping V^^; if U is smooth. 
We know that U is smooth by hypothesis that is normal concave utility. Therefore, for a 
given X* = \/U , the value U{x) — xy attains its supremum at U{x*). Therefore, the formula 
dMI is equal ioU{x*). □ 

Let [/ be a normal concave utility in T and a G M. Then, the continuous affine function 
X — > a — xy dominates U{x) if and only if 

Vx E T, a>U{x) — xy 

a > U{y) 
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Hence, we proceed to define the conjugate (or polar) of U by 

U{y) = sup{U{x) - xy}, y G M+ (45) 

a;>0 

which can be seen as the Lagrangian associated with the system of perturbations y £ M^. 

According to Rockafellar, in order that U is proper and U = U, it is necessary and sufficient 
that U is lower-semi continuous with respect to the weak topology induced on C by D, but 
we have already established continuity through the assumption that U belongs to the class of 
normal concave utilities. 
From (j45p . we obtain 

U{x) = inf [;7(y) +xy], x G M+ (46) 
From eqn. (j451 146p and lemma (j33p . we obtain the following inequalities 

U{I{y))>U{x) + y[I{y)-x],yx>0,y>0 (47) 
U{U'ix)) < Uiy)-x[U'{x)-y],yx > 0, y > (48) 
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B Static Optimization problem 

The set of wealth X^'^[-) financed by tt and having initial endowment x is given by 

V = {Tr€A: (X^'^) G Cl{P)} 

where £^(P) = C\{Vl,T,P) and < p < oo. 

Now, consider the following static variational problem 

subject to 

E[HtA\ = ^ 

In the following proposition, we show that by solving the static problem, we obtain the 
solution to the dynamic problem (|26|) . and we prove equivalence by contradiction. The argu- 
ment is that an agent prefers a vr solution to the dynamic problem, whose budget constraint 
is satisfied. But if the solution is budget feasible, then it will coincide to the solution of the 
static problem. 

Proposition 34. Consider ^ is a solution to the static variational problem and takes values 
in V. Then ^ is the solution to the dynamic problem. The converse also holds. 

Proof. Assume there exists a solution ^ financed by vr and initial capital x such that the 
following inequality holds 

E[U{i)] > E[Um 

But from (|24p. we have 

E[HtJ] < x 

and deduce that ^ is budget feasible, which contradicts the hypothesis that ^ is a solution to 
the static problem. The second statement follows from elementary arguments. □ 
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C Misc. proofs 



Proof, (lemma [T^ 

The proof will follow in two stages. First, we prove that Zt^ is a martingale/supermartingale 
and that E[Zt^] = 1, and then we show that Pn « (absolutely continuous). 

• By Taylor expansion applied on eq. (jlip and discarding the higher order terms of At, 
we get 

n-1 

Zt^ = l-Y,Zt,eiAWt, (49) 

k=0 

which is a local martingale. Let an increasing sequence of Markov times r„, such as 
the sequence {Z {tn l\Tn) , J^t„) is uniformly integrable, and Z{tn f\Tn) satisfies Zt^_^ = 
E[Zt„\J-t„_-^]- Then, from monotone convergence and Fatou's lemma, we obtain that 
Zt„ is a supermartingale. The assertion that E[Zt^]=l, follows from eq. (j49ll . subject 
to the initial condition Zt^ = 1. 

By application of Taylor theorem on A(Wt„Zt^), we have 

MWt^ZtJ = Wt^AZt,, + Zt^AWt^ + A[Wt„,AZtJ 

= -WtAOtAWt^ + Zt^AWt^ + Zt„p0t„At - Zt^pOt^At 
= -(dlWt^ + \)Zt^AWt^ 
which is a local martingale. Therefore, 

which proves that Wt^ has the martingale property. 

• Subject to the assumption Z{tQ) = 1 i.e. Zt„ = 1, we have that for any measurable set 
A G J^t„ ! we have from Radon-Nikodym that 

f,w{A) = Pniu : R e A) = [ Zt„iLo)dPn= [ EiZMlJ'l^^jdPn (50) 

Let the random process r]t^{Lo) = E[Zt^{Lo)\J-^] satisfy the Kolmogorov criterion 

E[\7jt^^At-r]tX]<C\At\'+' 
for some given constants a > and e > 0. 

Then, 'qt^[u)) has a modification ^'(7^4^(0;)), that is progressively measurable with respect 
to ^f^, which we denote by 

^{7^tA^)) = E[Zau)\T^J 
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and we rewrite the last part of eq. (j50p as follows 



■.ReA J A 

and by ([50]) . we have 



A 

Hence, we obtain the absolute continuity of ^-nd 

— = 



□ 



Proof, (theorem lisp We follow similar arguments to (j23l ). 
From eq. (fTTj) and (fT2]) it follows 



n—l , n—l 



Zt„ = exp{- J] eiAWt, + IY1 ^IP^t.At} (51) 

fc=0 fc=0 

By Taylor expansion theorem, and by discarding the higher order terms of At, we get 

a(^) =Ct.C^^t. (52) 

Making use of the assumption that we can replace the quadratic variation of AW by its 
expectation, we have 

n-l 



k=0 



From the assumption above and E[\\Pt^ — < oo, it follows that A y-^j permits the 
following representation (see theorem 5.13 in (1241)) 



„ 1 _ 1 ""^ _ ^ 

E[^\Fl!] = E[—] + Y.Ei^MFu]^Wt, 



^ n—l 
k=0 



□ 
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Proof, (theorem [21]) . The proof follows similar arguments to (jisl ) , and is restated here only 
for continuity reasons. 

From the budget constraint inequality (j24|) and (|47p . we have 

u{o>u{i) + y{x)Htjc-i) 
>-\uii)\-y{x)Ht^ 

The boundness condition follows from the monotonicity of 7^^^ and the supermartingale 
property of Ztj^ . 

Moreover, we have 

uio > u{x!f) + y{x)HtA^ - xff) 

Therefore, from the above inequality, we obtain the required result 

E[Um > V{x) 
The same result can also be obtained from inequality 

U{x) <U{y) + xy, \/x,y > Q (53) 

which is derived by rearranging (j45|) . Consequently, from (j24|) and the inequality (|53|) . we 
have 

C/(xf/) < U{y{x)Ht,) + j;(x)/7t^xf/ (54) 

Therefore, 

V{x)<V{y)+y{x)E[H,^xl^] 
< V{y) + y{x)x 

The above inequality turns to equality, if and only if (j28p holds. The optimality of vr'^ is 
implied by the existence of optimal solution for tt'^. 

□ 



26 



References 



[1] A. Bensoussan. Stochastic control of partially observable systems. Cambridge University 
Press New York, NY, USA, 1992. 

[2] J. M. Bismut. Growth and optimal intertemporal allocation of risks. Journal of Economic 
Theory, 10:239287, 1975. 

[3] J. C. Cox and C. Huang. A Variational Problem Arising in Financial Economics. Sloan 
School of Management, Massachusetts Institute of Technology, 1988. 

[4] M. G. Crandall and P. L. Lions. Viscosity solutions of hamilton-jacobi equations. Trans- 
actions of the American Mathematical Society, 277(l):l-42, 1983. 

[5] D. Cuoco. Optimal consumption and equilibrium prices with portfolio constraints and 
stochastic income. Journal of Economic Theory, 72(l):33-73, 1997. 

[6] J. Cvitanic and I. Karatzas. Convex duality in constrained portfolio optimization. The 
Annals of Applied Probability, 2(4):767-818, 1992. 

[7] J. Cvitanic and I. Karatzas. Hedging and portfolio optimization under transaction costs: 
a martingale approach. Mathematical Finance, 6(2):133165, 1996. 

[8] J. Cvitanic, H. Pham, and N. Touzi. A closed-form solution to the problem of super- 
replication under transaction costs. Finance and Stochastics, 3(l):35-54, 1999. 

[9] M. H. A. Davis. Lectures on Stochastic Control and Nonlinear Filtering. Springer- Verlag 
New York, 1984. 

[10] M. H. A. Davis and AR Norman. Portfolio selection with transaction costs. Mathematics 
of Operations Research, 15(4):676-713, 1990. 

[11] I. Ekeland and R. Temam. Convex analysis and variational problems [russian translation]. 
Mir, Moscow, 1979. 

[12] J. M. Harrison and S. R. Pliska. Martingales and stochastic integrals in the theory of 
continous trading. Stochastic Processes and their Applications, 11:215-260, 1981. 

[13] Y. M. Kabanov and G. Last. Hedging under transaction costs in currency markets: a 
continuous-time model. Mathematical Finance, 12(l):63-70, 2002. 

[14] G. Kallianpur. Stochastic filtering theory. Applications of Mathematics, 13, 1980. 

[15] I. Karatzas, J. P. Lehoczky, and S. E. Shreve. Optimal portfolio and consumption deci- 
sions for a small investor on a finite horizon. SIAM Journal on Control and Optimization, 
25(6):1557-1586, 1987. 

[16] L Karatzas, J. P. Lehoczky, S. E. Shreve, and G. L. Xu. Martingale and duality meth- 
ods for utility maximization in an incomplete market. SIAM Journal on Control and 
Optimization, 29(3):702-730, 1991. 

[17] I. Karatzas and S. E. Shreve. Brownian motion and stochastic calculus. Springer- Verlag 
New York, 1991. 



27 



[18] I. Karatzas and S. E. Shreve. Methods of Mathematical Finance. Springer, 1998. 

[19] H. Konno and A. Wijayanayake. Portfolio optimization under dc transaction costs and 
minimal transaction unit constraints. Journal of Global Optimization, 22(1):137-154, 
2002. 

[20] D. Kramkov and W. Schachermayer. The asymptotic elasticity of utility functions and 
optimal investment in incomplete markets. The Annals of Applied Probability, 9(3):904- 
950, 1999. 

[21] V. Krishnan. Nonlinear filtering and smoothing. Wiley New York, 1984. 

[22] H. Kunita. Stochastic Flows and Stochastic Differential Equations. Cambridge University 
Press, 1990. 

[23] P. Lakner. Optimal trading strategy for an investor: the case of partial information. 
Stochastic Processes and their Applications, 76(l):77-97, 1998. 

[24] R. Lipster and A. Shiryaev. Statistics of Random Processes. Springer, 2000. 

[25] H. Markowitz. Portfolio selection. The Journal of Finance, 7(1):77-91, 1952. 

[26] R. C. Merton. Lifetime portfolio selection under uncertainty: The continuous-time case. 
The review of economics and statistics, 51(3):247-257, 1969. 

[27] R. C. Merton. Optimum Consumption and Portfolio Rules in a Continuous-time Model. 
MIT, 1970. 

[28] G. Montana, K. Triantafyllopoulos, and T. Tsagaris. Flexible least squares for temporal 
data mining and statistical arbitrage, eprint arXiv: 0709.3884, 2007. 

[29] B. Oksendal. Stochastic Differential Equations: An Introduction with Applications. 
Springer, 2006. 

[30] H. Pham and M. C. Qucncz. Optimal portfolio in partially observed stochastic volatility 
models. The Annals of Applied Probability, ll(l):210-238, 2001. 

[31] R. T. Rockafellar. Integrals which are convex functionals. Pacific J. Math, 24:525-339, 
1968. 

[32] R. T. Rockafellar. Convex Analysis. Princeton University Press, 1996. 
[33] S. Shreve. Stochastic Calculus for Finance II. Springer, 2004. 

[34] S.E. Shreve and H.M. Soner. Optimal investment and consumption with transaction 
costs. The Annals of Applied Probability, 4(3):609-692, 1994. 

[35] G. L. Xu and S. E. Shreve. A duality method for optimal consumption and investment 
under short-selling prohibition, i. general market coefficients. The Annals of Applied 
Probability, 2(1):87-112, 1992. 

[36] Thalcia Zariphopoulou. Investment-consumption models with transaction fees and 
markov-chain parameters, 1992. 



28 



